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PROVISIONAL SPECIFICATION 

Invention Title: 

Quantum device 
The invention is described in the following statement: 
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Quantum Device 

Technical Field 

This invention concerns a quantum device* suitable for quantum computing, based on 
3 dopant atoms located in a solid semiconductor or insulator substrate. In further aspects 
the device is scaled up. The invention also concerns methods of reading out from the 
devices, initializing them, using them to perform logic operations and making them. 

Background Art 

10 There are a number of existing proposals for quantum computing (QC) based on charge 
qubits in the solid-state. There is also the Kane proposal for a real spin quantum device 
based upon phosphorus donor atoms in silicon. Proposals relevant to this invention are 
now presented in chronological order: 

1 S Logical devices Implemented using quanUmt cellular automata^ 
P,D. Tougaw and C.S. Lent, J. Appl. Phys. 75, 1818 (1994). 

Conditional quantum dynamics and logic gafes^ 

A. Barenco, D. Deutsch, Ekert and Josza, Phys, Re\K Lett. 74, 4083 (1995). 

20 

Minimal energy requirements in communication^ 
R. Landauer, Science 272, 1914 (1996). 

Coupled qumUum dots as quantum XOR gate^ 
25 J, A. Brum and P. Hawiylak, Superlattices midMicrostr. 22, 43 1 (1997). 

Quantum information in semiconductors: Noiseless encoding in a quantumniot array, 
P. Zanardi and F. Rossi, Phys. Rev. Lett. 81, 4752 (1998). 

30 

Quantum computation with quanttm dots^ 

D. Loss and D.P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 

35 A silicon-based nuclear spin quantuiu computer, 

B. E. Kane, Nature 393, 133 (1998). 
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Microwa\fe spectroscopy of a qiiautum-dot molecttle^ 

T.H. Oosierkamp, T. Fujisawa, W.G. van der Wiel, K. Ishibashi, R. V. Hijman, S. 
Tarucha and L.P. Kouwenhoven, iVa/?/re 395, 873 (1998), 

5 Dephasmg in open qucmUiw dots, 

A.G. Huibers, M. Switkes, CM. Marcus, K, Campman and A.C. Gossard, Phys, Rev. 
Lett 81, 200 (1998). 

Coupled quanhim dots as quanium gates^ 
1 0 G. Burkard, D. Loss and D.P. DiVincenzo, Phys, Rev. B 59, 2070 (1999). 



Oucaitum gates by cottpled asymwe/ric quatUum dots and CNOT-gale operation^ 
T. Tanamoto, Phys. Rev. A 61, 022305 (2000). 

15 

Coherent properties of a iwo^level systent based on a quantum-dot photodiode, 

A. Zrenner, E. Beham, S. Stufler, F. Findeis, M. Bichler and G. Abstrieter, NaUtre 418, 

612(2002). 

20 

High coherent solid-state qubit from a pair of quantum dots, 
Xin-Qi U and YiJing Yan, quant-ph/0204027 (2002). 

Two-electron quantum dots as scalable qubits, 
25 J.H. Jefferson. M. Fearn, D.L.J. Tipton and T.P. Spiller, quant-ph/0206075 (2002). 

The potential of quantum information processing relies on one of the most fragile 
resources in nature - quantum entanglement. Even without the destructive influence of 
environmental decoherence a practical quantum computer requires hundreds, if not 
30 thousands of gates to perform useful tasks. Since a system of qubits can never be totally 
isolated from the environment, error-correcting schemes must be incorporated in the 
design and encoding of the qubits, further increasing the number of qubits required. 

The real challenge is to prove the principle of QC in a technology which is inherently 
35 scalable. This does not simply mean the ability to make a lai^e number of 
interconnected qubits. The concept of scalability requires that the error probability per 
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qubit be kept below some threshold determined by fault tolerant implementations. 
What that threshold actually is depends on the particular implementation and the 
physical origin of decoherence. 

S Disclosure of the Invention 

The invention is a quantum device comprising a pair of dopant atoms located in an 
otherwise electrically inert solid substrate, such as a semiconductor or insulator, in 
which it is possible to ionise the dopant atoms. The dopant may be an electron donor 
10 or an electron acceptor. The two dopants produce a double-well electric potential and a 
charge qubit is realised by the location of one or more electrons or holes within this 
potential. 

In one example a pair of phosphorus (P) atoms are located in a silicon (Si) substrate. 

1 5 lonisation of one of the donor electrons creates a P-P"** system, equivalent to an artificial 
molecule, and the charge qubit is formed by the location of a single electron in one 
or other well of the double-well potential. The spacing between the P atoms may be up 
to 200nm. More specifically it may be in the range 20 to lOOnm. The basic qubit and 
definition of logical states can be expressed in terms of left and right localized states 

20 with respect to the double-well potential. 

The two P atoms may be buried up to 200nm below the surface of the device. More 
specifically they may be buried in the range S to SOnm. 

25 The silicon substrate may be coated with an insulating layer, such as SiOz, to isolate the 
donor electrons from any surface electrodes. 

Gate electrodes may be placed on the surface of the substrate above the donor atoms to 
allow for external control of the charge waveftinctions. A first gate may be located 

30 over the potential barrier, between the two wells, to control the barrier height, this gate 
is known as the 'barrier' gate or B-gate. Another gate, or more than one gate, may be 
provided for elsewhere to control the relative shapes and sizes of the two wells, this 
gate is known as the 'potential off-set', 'symmetry* or S-gate. Suitable biasing of these 
gate electrodes may allow one qubit logic operations, such as the quantum NOT 

3S operation, to be performed. 
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One or more charge detection devices may also be provided on the surface of the 
substrate for qubit readout and to confirm initialization of qubits. The charge detection 
device may be a sensitive field-effect transistor, or it may be a single electron transistor 
(SET). 

5 

The configuration described, comprising two P atoms, surface gates and read-out SETs 
may be repeated many times on a single silicon chip, allowing for scale*up to a many 
qubit processor. 

10 The entire device may be cooled well-below room temperature to ensure that the 
donors are not thermally ionised and to minimise coupling of the charge qubits to the 
environment. More specifically it may be cooled to liquid helium temperatures of 4K 
or below. 

15 In contrast with real spin qubits, a charge qubit can couple strongly to fluctuating local 
charge environments. This kind of coupling leads to fluctuations in the relative energies 
of the charge qubits and thus will be largely responsible for dephasing errors. 

In comparison with existing charge qubit proposals based upon quantum dots (listed 
20 above) the P-P"** charge qubit is the lowest filling quantum dot possible, with the 
simplest quantum-substructure. In particular, the potential wells that confine the 
electrons are created by single ionic potentials, rather than by lithographically patterned 
electrodes or clusters of atoms. For these reasons we expect the P-P'*' system to suffer 
less fi-om decoherence than the larger quantum dot charge qubits. 

25 

It is interesting to compare directly the real-spin Kane proposal [1] with this pseudo- 
spin P-P* proposal, [n the charge qubit case we sensibly want the system to be cooled 
as far as is possible to minimize environmental coupling, but we do not necessarily 
require a magnetic field. The decoherence of the charge qubit system will be much 
30 faster than the real-spin qubit, however, the intrinsic "clock speed" of the charge qubit 
is also much faster. In the case of the Kane design the CNOT gate is implemented 
adiabatically - the clock speed is therefore restricted to the nuclear Rabi fi'equency, of 
order kHz. In the P-P* case the barrier gate may be used to tune the natural clock speed 
down to the level of fast pulse technology, of order GHz. 

35 



Since we are no longer dealing with nuclear spin states, isotopic impurities (e.g. Si^^) 
will not aflfect the charge qubit. Therefore, one can use standard Si wafers in the 
fabrication process. Although we will use the P-doped Si model throughout this 
proposal, there is no intrinsic restriction to the P:Si system. Alternatives such as Si- 
S doped GaAs would also be possible. 

Furthermore, since the Coulomb interaction acts over a longer range than spin-spin 
interactions, the charge qubit devices can be made larger than their spin-based 
counterparts. In particular, P donors in the Kane architecture are ideally placed 20nm 
10 apart, while the donor separation for charge qubits could be as large as 200nm with the 
appropriate choice of substrate and gate potentials. This increased scale makes 
nanoiabrication considerably easier. 

Given that the fabrication process is very similar, the main advantage of the real-spin 
IS over the charge qubit is the long decoherence times. However, the trade-ofF is in the 
measurement process: measurement of the charge state qubit is considerably easier to 
perform than the single-spin measurement scheme as envisaged by Kane. 

20 In another aspect, the invention is a quantum device comprising two pairs of dopant 
atoms, each located in a solid substrate such that two double-well potentials are created, 
and two charge qubits are formed by the location of one or more electrons within the 
double-well potentials. And where, one or more gates above the substrate control 
electron tunneling between the two qubits. And where, by suitable biasing of these and 

25 other gate electrodes, two-qubit logic operations, such as the quantum SWAP 
operation^ may be performed. 

In another aspect, the invention is a quantum device comprising two pairs of dopant 
atoms, each located in a solid substrate such that two double-well potentials are created, 
30 and two charge qubits are formed by the location of one or more electrons within the 
double-well potentials. And where, one or more gates above the substrate controls the 
Coulomb interaction between the two qubits. And where, by suitable biasing of these 
and other gate electrodes, two-qubit logic operations, such as the quantum CNOT or 
CPHASE operation, may be performed. 

35 
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In another aspect, the invention is a quantum device comprising more than one of the 
devices comprising two pairs of donor atoms, arranged adjacent each other, and having 
gates to enable interaction between each adjacent pair. 

5 In a further aspect the invention is a method of reading out the logic state of a qubit 
from such a quantum device, comprising the steps of: 

(i) Biasing the S-gate to a non-zero value to define a known relaxed state; 

(ii) Turning on the SET; 

(iii) Waiting for a period of time for the system to relax to the relaxed state; 
10 (iv) Monitoring SET current during the waiting period to determine whether 

a change of state occurs, or not; 
(v) Inferring the state before step (i) from the behaviour of the SET current 
monitored during step (iv), 

IS In a further aspect the invention is a method of initializating the logic state of each 
qubit in such a quantum device, comprising the steps of: 

(i) Biasing the S-gate to a non-zero value to define a known relaxed state; 

(ii) Waiting for a period of time for the system to relax to the relaxed state. 

20 In a fiirther aspect the invention is a method of using such a quantum device to perform 
quantum computations, comprising typical steps which include: 
CO Initializing each qubit in the quantum device; 

(ii) Using a suitable quantum algorithm to carry out one-qubit and two-qubit 
operations on the qubits in the device via an appropriately timed 

25 sequence of control gate biasing; 

(iii) Reading out the logical states of the qubits. 



In a further aspect the invention is a method of making such a quantum device, 
comprising the steps of: 
30 (i) creating a resist mask for implantation of an array of one or more pairs 

of dopant atoms into a semiconductor substrate; 

(ii) implanting one dopant ion through each hole in the mask to form an 
array of buried-atom qubits; 

(iii) fabricating conducting electrodes above each implanted pair of donor 
35 atoms to serve as a B gate, an S gate and a SET; 
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(iv) annealing the device to a temperature sufficient to remove implant 
lattice damage and activate the dopant - this step may occur either 
immediately before or after step (iii); 

(v) pre-^initializing each pair of donor atoms to remove one of the electrons. 

5 

In a further aspect the invention is a method of making such a quantum device, 
comprising the steps of: 

(i) passivating a clean semiconductor surface; 

(ii) patterning the passivated surface to open an array of holes; 

10 (iii) exposing the resulting surface to donor atoms, perhaps in the form of 

dopant-bearing molecules; 

(iv) overgrowing the dopant atom array with further layers of the 
semiconductor and annealing if necessary; 

(v) coating the surface with an insulating layer if necessary; 

IS (vi) fabricating conducting electrodes above each implanted pair of donor 

atoms to serve as a B gate, an S gate and a SET; 
(vii) pre-initialLzing each pair of donor atoms to remove one of the electrons. 



Brief Description of the Drawings 
20 Some examples of the invention will now be described with reference to the 
accompanying drawings, in which: 

Figure 1(a) is a pictorial view of aP-P* charge qubit; Figure 1(b) is an elevational view 
of the charge qubit; and Figure 1(c) is a plan view of the charge qubit. 

25 

Figure 2(a) is a potential diagram showing the double-well potential created by the 
charge qubit; Figure 2(b) is a potential diagram of the P-P* charge qubit with a black 
dot representing a localised electron state in the left hand well of the double-well 
potential; and Figure 2(c) is a potential diagram of the P-P^ charge qubit with a black 
30 dot representings a localized electron state in the right hand well of the double-well 
potential. 

Figure 3(a) is a potential diagram of the P-P^ charge qubit showing the effect on the 
barrier of a non zero voltage on the B*gate; Figure 3(b) is a potential diagram of the P- 
35 P* charge qubit showing the effect on the double-well of a non zero voltage on the S- 
gate. 
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Figure 4(a) is a potential diagram of a symmetric molecular orbital wave function; 
Figure 4(b) is a potential diagram of an anti-symmetric molecular orbital wave 
function. 

5 

Figure 5 is a graph showing the ungated tunnelling rate, expressed as frequency f, 
against the separation R, in units of Bohr radius, a« 

Figure 6 is a graph of the wave functions energies (E) against bias on the S-gate, and 
1 0 includes associated diagrams of the effect on the double-well potential. 

Figure 7(a) is a logarithmic plot of the amount of non-logic state contamination in the 
ground and first excited states for a donor spacing R=6.5a; Figure 7(b) is a logarithmic 
plot of the amount of non-logic state contamination in the ground and first excited 
1 5 states for a donor spacing R=9a. . 

Figure 8 is a graph showing Energy gap vs. AB around point ^o~0,46884, 5b=0. 

Figure 9 is a plan view of coupled qubits in the SWAP configuration. 

20 

Figure 10(a) is a plan view of coupled qubits in the CNOT configuration; Figure 10(b) 
is a diagram showing the effect of the state of qubit Q, on the barrier height of Q2 

Figure 1 1(a) is a plan view of coupled qubits in the staggered configuration; qubit Qi 
25 on the potential symmetry seen by Q2. 

Figure 12(a) is a two dimensional representation of two coupled ungated qubits in an 
arbitrary configuration; Figure 12(b) is a two dimensional representation of two 
coupled ungated qubits in the CNOT configuration; and Figure 12(c) is a two 
30 dimensional representation of two coupled ungated qubits in the CPHASE 
configuration. 

Figure 13 is a plan view of N-qubit architecture based on the SWAP configuration. 

35 Figure 14 is a plan view of N-qubit architecture based on the staggered CPHASE 
configuration. 
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Figure 15 a plan view of one of many possible N-node quantun) cellular automata 
networks, which could be constructed using the P-P^ qubit. 

5 Figure 16 is a schematic flow chart showing the readout/initialization process. 

Figure 17(a) is a diagram showing the effect of a varying S-gate potential on the 
eigenstates of a single qubit; Figure 17(b) is a pulse-timing diagram showing the 
sequence of pulses on the S-gate and the SET in order to study coherent one-qubit 
10 rotations. 

Figure 18 is an SEM image showing two SETs (a "twin-SET' device) which can be 
operated at a frequency up to 500MHz, which is the frequency range required to 
perform charge qubit read-out. 

IS 

Figure 19 is experimental data obtained with a device like that shown in Figure 18. 
The data shows the current through SET-1 and SET-2 as single electrons are transferred 
between two metallic islands located between the SETs. 

20 Figure 20 is experimental data obtained with a device like that shown in Figure 18. 
The SET can respond to a changing potential (or charge movement) on a timescale of 
100ns» which is expected to be less than the relaxation time Ti for the qubit. This 
indicates that charge qubit read-out should be possible with such a device. 

25 Figure 21(a) is a schematic of a device with two implanted donors> two read-out SETs, 
and a series of five gates which can be used as B-gates and S-gates to control the 
donors; Figure 21(b) is an SEM image of such a device that has been fabricated using a 
top-down production process, in which clusters of P ions have been implanted at the 
sites marked with circles. This device can be used for charge qubit proof-of-principle 

30 experiments. 

Best Modes of the Invention 
The P-P*^ charge qubil 

35 
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Referring first to Figures I and 2, a charge qubit 10 is based on electrons localized by 
two donor P atoms 11 and 12 buried in a silicon substrate 13. Each qubit comprises an 
"artificial molecule" P-P"^ system in which a single electron 14 remains, the other donor 
electron must previously be removed from the system by an ionization process, such as 
5 electric filed induced ionization produced by suitable gate biaising. Gates placed above 
the buried P-P* system, on a non-conducting barrier layer 1 5 allow for external control. 

The remaining electron 14 exists in a double-well potential 16, in one or other of the 
wells. In the potential diagrams the double-well potential 16 has a left-hand well 17 

10 and a right hand well 18. Between the wells is a potential barrier 19, and there are 
potential barriers 20 and 21 outside the wells. A basic qubit and definition of logical 
states in terms of left 17 and right 18 localized states with respect to the double-well 
potential 16 is shown in Figure 2(b) and Figure 2(c), respectively. In order to facilitate 
discussion of coupled qubit systems we also introduce the occupation number 

1 5 definition such as ID = Lj 1 0> = 1L,0> etc. 

A barrier gate (B-gate) 22 placed between the two P atoms is used to exert control over 
the height of barrier 19. A potential off-set, or symmetry (S-gate) 23 is placed nearby . 
A SET 24 is provided for readout. 

20 

Figure 3 illustrates the basic external control over the intra-qubit potential via the B and 
S gates. Figure 3(a) shows the effect of a non-zero voltage on the B gate, note in 
particular that barrier 19 is raised. Figure 3(b) shows the effect of a non-zero voltage 
on the S gate, note in particular the raised right-hand shoulder 21, the asymmetrical 
25 barrier 19 and the reduced left-hand shoulder 20. 

The P-P* charge qubit is the lowest filling quantum dot possible, with the simplest 
quantum-substructure. For this reason we expect the P-P^ system to suffer less fi-om 
decoherence than the larger quantum dot qubits. On the other hand, a charge qubit can 
30 couple strongly to fluctuating local charge environments. This kind of coupling leads to 
fluctuations in the relative energies of the charge qubits and thus will be responsible for 
dephasing errors. 

In the pseudo-spin P-P*** proposal, the system should be cooled as far as is possible to 
35 minimize environmental coupling. The decoherence of the charge qubit system will be 
fast, however, the intrinsic "clock speed" of the charge qubit is much faster. The barrier 
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gate will be required to tune the natural clock speed down to the level of fast pulse 
technology, of order GHz. 



Since we are not* dealing with nuclear spin states, isotropic impurities (e.g. Si^^ will not 



S affect the charge qubit. Therefore, one can use standard wafers in the fabrication 
process. 

Single Qubit Dynamics 
10 Ungated qubit 

In this section we discuss control over single qubit dynamics via external gate biases. 
Essentially we need to be able to initialize the qubit and perform arbitrary rotations. 
The first problem is the fast natural tunneling rate of the system for typical donor 
1 S separations. This is also a problem for two-qubit gates. 

To get first order estimates of the tunneling rates and gate control over the dynamics for 
the single qubit system we consider the P-P^ system using the usual effective Bohr 
radius of = 3nm appropriate to silicon. In the case of zero gate potentials the 
20 Hamiltonian in this effective-mass approach is taken to be: 



where the origin is placed mid-way between the donor nuclei along the x-^axis. The 
2S eigenstates of this Hamiltonian are known exactly, so we have constructed the basis of 
states {(Piif \R) } with i=l.J/ (the dependence on separation distance is manifest). In 
Figure 4(a) and (b) shows the first two lowest energy states, ^s(^^^) 9^4 ('^i^)! 
corresponding to symmetric and anti-symmetric configurations respectively. In the 
absence of external gate potentials the energy difference between these states 
30 AE{R)^EsiR) - Ea^H) determines the tunneling rate. The dependence on the tunneling 
rate with donor separation for zero gate bias is shown in Figure S - normally this rate is 
very fast, but with appropriate gate bias this system can be controlled. 



^0 




|r-7?/2| lr+/?/2|' 



In terms of these two lowest zero-gate eigenstates the computational space is given by: 

35 
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At zero-gating, the logic states are given in terms of the eigenstates of the system as: 

5 

|o> <^ vl in - -j^iiPs(r)^<PA(f)), 

Gated Qubit 

10 

Now we consider the placement of gates according to the single-gate configuration of 
Figure 1(c): a B-gate 22 controlling the barrier height, and as S-gate 23 which allows us 
to raise and lower the potential asymmetrically. The spatial dependence of the gate 
potentials, and Vgir)^ in the silicon substrate can be conveniently calculated 

15 using the TCAD Poisson solver (gate geometry, widths etc). In effect we add to the 
Hamiltonian terms of the form: 

20 where the time-dependent pulse profiles SQ) and B(t) describe the time evolution 
(scaling) of the gate biases. 

The struaure of the lowest two states in BS space is complicated. Around the zero-gate 
point, i5— iS-0, the situation is fairly clear, and is shown in Figure 6 as the conditions for 
25 minimum energy gap between the first two states. The ground state is this symmetric 
state and the excited state is the anti-symmetric state. As we change the S-gate voltage 
adiabatically the states evolve to localized states as shown to the right and left. 

Loss in fidelity due to ''state-leakage" will occur if the logic states (localized states) are 
30 not exact eigenstates of the Hamiltonian. This effect should be minimized by 
positioning the P-donors far enough apart, and applying the appropriate S-gate vohage. 
In the current qubit model, the variation of the error due to state-leakage is plotted in 
Figure 7 for a range of S-gate voltages. Clearly, for separations of 9a or greater. Figure 
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7(b), this source of fidelity loss is less than the 10"** benchmark. At So = 0 the first two 
states are equal combinations of the localized states. At |So| = 0,3 volts or higher the 
qubits are quite accurately defined by the localized states. 

5 In order to be able to describe the dynamics of the gated P-P*** qubit we need the 
effective Hamiltonian in the computation sub-space. First, we consider the eigenstates 
of the system when the gate strengths B and S are non zero. From the exact treatment of 
the unperturbed system we construct the Hamiltonian matrix with respect to N basis 
states to obtain the matrix: 

10 

^ejriR\j ^<<PiinR)\H\pj(f;R)>. 

Diagonalization of this matrix for non-zero gate strengths gives a new energy gap 
between the lowest states AE{R,S,B). This is of course a non-trivial functional 
15 dependence given that the potentials are described by TCAD simulation in the first 
place. Nevertheless, we have solved the system for R-6,5a and 7?=9.0a for various gate 
potentials. 

Effective Hamiltonians and Qubit Rotations 

20 

\Ve need to find the regions in (B,S) space where single qubit operations in the lowest 
two states are feasible, that is the tunneling rate can be controlled and slowed down. To 
this end we compute the effective Hamiltonian for deviations (^(0. ^(jO) about a 
given point {Bo^So). In general this will be of the form: 

25 

where 

K (0 = (So , fio ) + .^s (So , Bo )m) + (^0 . Bo ) AB(0 . 
hAO^ZoiSo.Bo)-^Z^iSo.Bo)ASii) + Z^(So.Bo)mO 

30 

The all-pervading dependence on R has been suppressed. In Table 1 we list the vahies 
of the effective Hamiltonian coefficients calculated in the hydrogenic approximation. 
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This form for Hq makes it convenient to study the dynamics of the qubit using NMR 
type pulses, where the rise time of the pulses are very short in comparison to their 
duration. Clearly, we can see how the function hj^i) controls qubit rotations and its 
composition in terms of the S and B potentials. No matter how we define our 
5 computational basis, a n/2 rotation exp(-i;ro^2) is implemented by balancing the 
coefficient of oi to zero and then allowing the system to evolve for a time Top widi the 
(Tx coefficient non-zero. Changing the basis from localized states to delocalized states 
induces the transformation ok so care must be taken on this point. 

10 The main question is the limitation on the operation time - the natural frequency of the 
un-gated system (Figure 5) is very fast for donor separations up to R=10a so we need to 
find gating operations to slow the tunnelling rate down. We can achieve the balance 
conditions by suitable control over and AS(/) pulses: ^^(A5*,A5*)sOand 

^y(Ai5\A5*) sO. This gives: 

IS 

The values of AB* and AS* for three gate configurations {Bq,So) are given in Table 2. 
We examine first how a n/2 rotation exp(-/;rG5^/2) may be applied by B-gate control. 
20 This is carried out by maintaining hz=0 for AB(t)^ AB* using: 

^^ -(Zq^Z^AB) ^ 

which gives: 

hAm^~[{XoZs-ZoXsh{X^Zs+Z^Xs)ABl 
Zs 

25 The timing of the AB pulse can now be read off from: 

Top =^[A,(AB)J-' =^^Zs[{XoZs-ZoXs)+{X^Zs+Z^Xs)AB]'\ 

(recall that is in units of 45meV). In Table 3 the voltage pulse magnitudes required 
30 for the different gate configurations (Bo.So) and donor separations (R=6.5 a and R=9a) 
are given for what might be considered extremal values of pulse times: Top = 1 ns and 
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fop = I //s. In each of the cases the voltages required are feasible, demonstrating that the 
single qubit operation time can be controlled 

In Figure 8 we show the situation around the configuration 5(rM).46884, 5H>. 

5 

The non-zero B-gate case deserves some comment. First, in an ideally symmetric 
qubit/gate geometry, one does not expect a minimum to develop as the barrier height is 
increased when S=0. In fact as the barrier is raised, by making B(t) more negative, the 
response of the first two states will be different - the dependence of the energies of the 

10 first two states on B-gate voltage is to first order linear, but with different slopes. 
Hence, at some point the energy levels will cross. Any asymmetry in the system will 
effectively introduce a coupling between the states creating an avoided crossing. The 
calculations involve rounding errors, and so the cross-over becomes an avoided 
crossing, and we obtain the behaviour as seen in Figure 8. A straightforward LCAO 

IS calculation using a symmetric barrier potential confirms this. 

Indeed, in an actual device one does not expect to obtain the perfectly symmetric case, 
and a certain amount of device characterization will need to be carried out to determine 
the adjustment required on the gate potentials to facilitate operation, as per the 
20 estimates given above. 
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Table 1: Single qubit efTective Hamiltonian coefilcients for various gate 
conHgu rations 
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5b=0 
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-1.2228 
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-0.7714 
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? 
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^0 




-0.1417 




? 


0.4663 




? 


-0.0032 
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Table 2: Gate configurations and balance conditions for single qubit rotations. 



Gate confij^ration 


R=6.5a 


R=9.0a 


5o=0, 50=0 




Ai?' = -1.00585 




^5^ = 


AS' = -2.1x10-* 


-\/2 








Bo=-0.46884, Sfi=0 




AB' = +0.003683 












A6- = 


A^* =+0.000101 


B,f=0, 5b=-0.3 


A5' = 


AB* = -0.533922 










92 =«'R 




A5' = +0.300349 


Table 3: Single qubit rotations: pulse timing 


Gate configuration 


R=6.5a 


R=9.0a 




Top =0. Ins 




r<,^=0.1ns 








AS = -1.0060 


i5o=0, 5o=0 


AS = 




A5 = +0.0005 


92 = -j^{9L-9R) 


rq,=lns 




ro^=lns 






Ai5 = -1.00587 


A5 = 




A5 = +4.59x10-' 




Top =0. Ins 




ro^=0.1ns 




A5 = 




AB = +0.0336 


J?o=-0.46884, 5o=0 






AS = -0.000 1 


9i =9'l 








9^ =«»R 


r^p=lns 
A5 = 




r«p=lns 
A5= +0.00667 
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A5 = +7.66x10-* 




T„^=0.1ns 






AB = 


A5 = -0.4752 


5o=0, 5o=-0.3 


AS = 


AS = +0.3008 


















AB = -0.52805 






M= +0.30039 



Two Qubit Gates 



The coupling of qubits can proceed in a number of ways. Three basic coupling 
schemes, the first involving a direct external-gate controlled tunneling interaction - the 
5 "SWAP" configuration - and the other two an indirect qubit-controlled coupling based 
on the Coulomb interaction between neighbouring qubits - the "CNOT* and 
"CPHASE" configurations - are shown in Figures* 9, 10 and 1 1 respectively. 

SWAP Gate Operation 

10 Referring to Figure 9, two pairs of charge qubits 30, 31 in the SWAP configuration are 
located parallel to each other in a substrate. Each individual qubit 30, 31 is constructed 
as set out in the description of the P-P"*" charge qubit above, each having a single 
electron 32, 33 remaining in the P-P^ system, a B-gate 34, 35, a S-gate 36, 37 and a 
SET 38, 39. The diagram shows the S-gates 36, 37 and the SETs 38, 39 in diagonally 

1 S opposite ends, but they could equally both be on the same side. 

Inter-qubit gates are also provided to control electrons tunnelling between the two 
qubits. Two further B-gates 40, 41 are provided between and parallel to the two qubits 
30, 31. Together, these B-gates 40, 41 extend the length of the qubits 30, 31 and their 
20 associated gates, except for a via-gate 42 that separates them and lies directly between 
the two qubits 30, 31. The via-gate 42 prevents diagonal transitions between the two 
qubits 30, 31. 

These gates are relatively straightforward to understand at the level of the single qubit 
25 dynamics described above (modulo the double |UR> occupancy step). Essentially the 
swap occurs by controlling the inter-qubit barrier gates & (40) and B4 (41) while the 
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intra-qubit dynamics are static by balance of Bi (34)and B2 (35). Assuming we are 
worldng in a left-right localization basis (w.r.t. the x-axis) the qubit computational 
states are labelled (interchangeably): 

|1) H/?) = M>" 

The notation |L,0) and |0,R> pave the way for a doubly occupied qubit |L,R>, which 
will occur as an intermediate state in the SWAP process. Starting with an arbitrary state 
of the two qubit system, j^,) = ®|^«/^)2>we wish to effect the SWAP operation on 

10 its components: 

\op) ^ |0)|0) 

m m 
\m ^ |o)|i)- 



By appropriate gate conttx)!, the evolution of the initial the two qubit state, 
1 5 I %) = 1^), ® 1^)2 . occurs in the following sequence (schematic): 



a,a2\L,Q\\L,0)^ + a,/?2(z:,0),|o,i?)2 + fi^a2\0,R\\L,0)^ + p,fi2\0,R\\0,Ii)2 
Bji B^i Bji B^-l 

^ a,a2lL.0),|L.0)2 + aMJ^.R^lOfi)^ + M2\0,0\\L,R)2 + A^2|0.^>i|0.^>2 
B^i B^i ^4 4 54 4' 

-> a,tt2|L.0>,|£,0)2 4- aMO,R\\L,0)^ + /?,a2|Z..0>^|0./?)2 + Ai^alO.^), 1 0.-^)2 

l^/> = aia2|i.0),l40)2 +aMO,R\\L,0)^ +/?,a2|^.0),|0./?)2 + AA|0./?),|0.;?>2 



By symmetry one expects transitions such as 

20 

|o./e),|o.;e>^ lo.o>,lo.2i?)2 




not to occur: although the actual electrons in this state will be exchanged when B3 is 
switched on, the product state remains effectively the same. If there is some residual 
asymmetry which generates this transition, a background magnetic field could be used 
S to polarize the electrons and stop this transition through Pauli blocking. 

Coulomb-^oupled gate operations 

Figure 10 represents a configuration which is based on an indirect qubit-controlled 
coupling. In this CNOT configuration the two qubtts 30, 31 are oriented perpendicular 

10 in a substrate. Again, each individual qubit has an associated electron 32, 33, B-gate 
34, 35, S-gate 36, 37 and SET 38, 39 that perform the same functions as previously 
described. The SETs are shown in diagonally opposite comers of the configuration but 
this is not a requirement. The B-gate 34 of the first qubit 30, and the S-gate 37 of the 
second qubit 3 1 occupy the remaining diagonally opposite comers, with the S-gate 36 

1 5 of the first qubit 30 and the B-gate 35 of the second qubit 3 1 now in parallel. Note that 
the SET 39 and S-gate 37 of the right qubit could be swapped without loss of 
generality. 

A further B-gate 40 is provided to control electrons tunnelling between the two qubits. 
20 This B-gate 40 extends the length of two charge qubits 30, 31 and their associated 
gates, and lies parallel to the second qubit 3 1 . 

Figure 1 1 represents a fiirther indirect qubit-controlled coupling configuration* In this 
CPHASE configuration the two qubits 30, 3 1 and their respective gates are arranged as 
25 described in Figure 9, but the B*gate provided to control electrons tunnelling between 
the two qubits is configured as described in Figure 10 and the two qubits are ofiEset (or 
staggered) vertically with respect to each other. 

The CNOT and CPHASE gates, as defined above, rely on the Coulomb interaction 
30 between electrons in the qubits. Understanding the gate operation quantitatively 
requires a lengthy computation of the four-donor, two-electron system. However, we 
can make good estimates based on some simple constructions. 

First we take the qubit positions as defined in Figure 12 (x-y plane only for clarity) 
35 below: in the ith qubit, the jth donor atom is located at = (^ff^J'f/^^^) • The effective 

molecular Hamiltonian for the two electron system is then: 
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5 where 




10 

In order to obtain gate operation characteristics^ at least to the level of the effective 
hydrogenic q)proximation used in the single-gate case, we would need to diagonalize 
the Hamiltonian as given, with respect to a complete two-electron basis. An outline of 
IS these calculations is given below, however, we note here that such calculations, while a 
useful guide, are probably not going to give very much quantitative insight into the 
actual charge qubit coupling taking place in the solid-state environment. For that one 
needs the full Bloch-wave treatment, which is beyond the scope of this proposal. 

20 The simple LCAO approach is to defme a basis of states in which each electron is 
localized about one of the qubits (and is therefore distinguishable). Symmetric and anti- 
symmetric single electron qubit wavefunctions can be defined as: 

V2(1±A<'>) 

25 

where the overlap integral, A^'^ is given by (1 S orbitals real): 
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The basis for the two-qubit system { 4^^ : a==l .4} is then constructed from these single 
electron qubit states: 

5 

The two-electron molecular Hamiltonian, in the presence of B and S is then 
diagonalised with respect to this basis. Details of the gate operation, pulse timing and 
so forth can be computed in ever increasing levels of approximation. 

10 In lieu of this lengthy exercise, we ^ve an estimate of the coupling expected for 
various inter-qubit separations from a semi-classical calculation. In the localized basis 
a potential affecting the barrier height as a B-gate corresponds to cr, 
pseudo-spin operator, while S-gate type eflfects correspond to a^. Hence the effective 
interaction between qubits in the CNOT configuration of Figure 12(b), and the 

15 CPHASE configuration of Figure 12(c) will be dominated by /^NorOk^^^ <^^^ and 
/cphaseCBs^^^ interactions respectively. A preliminary calculation based on bringing a 
positive charge of \e to dn = 10-20a from the qubit ^ves a value of 10"^ - 10"^ On 
45 meV units). This in turn gives a coupled qubit operation time of T12 « 0.1-lns. 

20 Scaleable N-Qubit Architectures 

Based on the two gate configurations presented above, we generalize the SWAP and 
CPHASE gates to scaled up N-gate architectures. . 

25 N-SWAF gate architecture 

Figure 13 represents a N-qubit architecture, where the qubits 30 are in the N-SWAP 
configuration with each qubit 30 oriented parallel in a substrate to create a row of N 
qubits. Each individual qubit 30 is constructed substantially as set out in the 
description of the P-P* charge qubit above, with each qubit 30 having a single electron 

30 released from a P atom, a S-gate 36 and a SET 38. Each pair of S-gates 36 and SETs 
38 extend out of opposite ends of one qubit 30 in parallel with it. Between each pair of 
adjacent qubits, and on the ends of the row, a B-gate 34 is provided that extends the 
length of the qubits and their respective gates. A via-gate 42 is provided in the centre 
of the B-gate 34 and is in line with the centres of the qubits 30. A fiirther via-gate 43 is 

35 provided in the centre of each qubit 30. 
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N-CPHASE gate architecture 

Figure 14 also represents a N-qubit architecture, but the qubits 30 are in an alternate 
staggered CPHASE configuration. This configuration is substantially identical to that 
described in Figure 13, but with two differences. Firstly, the qubits 30, while still in 
5 parallel are now arranged on two different levels with the qubits 30 in the row now 
alternating between the two levels. Secondly, no via-gate is provided in the centre of 
the B-gates 34. 

Quantum Cellular Automata network 

ID An example of a N-node quantum cellular automa network is depicted in Figure 15. 
Each individual qubit 30 is constructed substantially as set out in the description of the 
P-P^ charge qubit above, v/ith the addition of a via-gate 43 provided in between the two 
P atoms of each qubit 30. The two rows of qubits 30 are placed on top of one another 
such that the qubits are also vertically aligned. In this formation, each qubit is placed 

15 in the substrate perpendicular to each qubit that it is adjacent to and is not diagonally 
adjacent. On either side of each qubit 30, a B-gate 34 is also provided. 

The inherent scalability of the charge qubit (and indeed the spin qubit) allows one to 
define a N-qubit quantum cellular automata (NQCA) network, as shown in Figure 15. 
20 The operation of such a device is not as clear as the NCPHASE or NSWAP networks, 
however, the design permits general computation in terms of generation of 
entanglement followed by carefully defined measurement sequences. NQCAs might 
also be useful in defining heuristic quantum algorithms in local analogies to classical 
neural networks. 

25 

Decoherence: dissipation and dephasing 

One of the advantages of spin qubits over the pseudo-spin charge qubits is the longer 
decoherence times expected of spin degrees of freedom in the solid-state. The charge 
30 qubit is expected to decohere much faster - however, as the analysis in this proposal 
have shown the inherent gate operation times are also commensurately faster. Below 
we give some estimates for times associated with various sources of decoherence, and 
show that the condition T^co))'^cp^^^ quantum computation can be feasibly maintained. 

This conclusion is supported by the measurement of the decoherence time in GaAs 
35 quantum dots of about 1 ns, which can be considered to be a strong lower bound on the 
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decoherence in the P-P^ case given the fact that the charge qubit is less e>q>05ed to 
environmental fluctuations than GaAs quantum dots. 

Phonon mechanisms: 

5 

For typical operational time scales (ns) the energy gap between the qubit states is about 
10"^ meV, which is comparable to the ambient thermal energy of k^T= 0.0086meV at T 
= lOOmK. However, only phonon modes with wavelengths of the order of the electron 
orbit will couple appreciably. For intra-qubit donor distances (defining the electron 
10 orbit size) of /? = 10a to 50a, the corresponding phonon energies are about 10"* times 
larger: E^onan = 13 - 67meV. The occupation of these modes at T = lOOmK will be 
negligible, and so presumably will be their effect on dissipation. 

An estimate of the effect of phonons on dephasing can be obtained based on the qubit 
IS interacting with an environment of acoustic phonons - represented by a bath of 
harmonic oscillators O.eggett et al, RNdP 1987 etc). We follow the analysis of 
Tanamoto (PRA 61, (2000) p02230S) for the case of quantum dots in a-Si02. The 
spectral function J(w) has contributions from Ohmic and super-Ohmic parts; the latter 
of which is probably the most relevant at T=100niK, and is given in the Debye 
20 approximation by: 

For a tunneling frequency of Ob an estimate for the fiistest damping rate of the under- 
damped coherent oscillations with the bath at T = 0 is: 

25 

where the parameters for Si are y = 9.5 eV (deformation potential)^ p - 2329 kg/m^ 
(density), and c = 9.04x10^ m/s (speed of sound). The shortest decoherence time (i,e. 
30 not using the renormalized tunneling rate for lack of microscopic physics) is estimated 
from Tphonon = ^f^so' Taking a bare tunneling time of T^p = 0.1 ns as our base 
operational time, we obtain Tphonon 324 ns, indicating that a large number of gate 
operations could be performed before phonon induced decoherence becomes an issue. 
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Gate fluctuations and Johnson noise: 

The effect of gate fluauations is expected to be another source of decoherence for the 
S charge qubit. An analysis of Johnson noise can be carried out using the Master equation 
approach [ref Wellard and HoUenberg]. Fluctuations in the S and B-gate voltages are 
introduced stochastically as: 



10 



IS 



20 



25 



30 



The factors Agand scale the fluctuations and are given by the Johnson formula 
(assuming effective resistance and temperature of the gates are the same): 

Taking a n/2 rotation (dominated by the S-gate voltage) as a typical example of a gate 
operation, the dominant off-diagonal contribution to the density matrix is 

from which we obtain an estimate of the dephasing time due to Johnson noise in the 
gates: 

2h^ n 4.3x10"^ 



'Johnson r v2 



Rk^r Xl RT 

It should be noted that typically, a full-scale calculation of the decoherence^ which does 
not rely on the Master equation approach would get a softer dependence for short time. 

Charge traps: 

As a result of the fabrication process the inadvertent presence of charge traps is another 
possible source of unwanted environmental interaction. In particular, traps at the Si/SiOa 
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interface may be problematic. In the top-down process, the implantation of the P donors 
will also create traps, however, annealing after P implantation is expected to reduce the 
density of crystal imperfections to a level where the mean distance between charge 
traps is much larger than the intra-qubit distance. In many cases, the presence of a 
5 charge-trap near the qubit will give rise to a loss of fidelity rather than decoherence^ 
since the trap is essentially a static modification of the crystal potential. Problems 
remain, however, if tunneling between traps, or between traps and gates or SETs, can 
occur, since this dynamic process will lead to decoherence. 

10 Readout and Initialization 

A basic process for both readout and initialization is outlined below and summarized in 
Figure 16. The S-gate 23 is biased to a non-zero value, to define a well known relaxed 
state 50. The SET 24 is then turned on 51, and the system will either fall in the lefk 17 

15 or right 18 of the localised states of the double-well potential. While waiting a certain 
period 52 for the system to relax to a relaxed state, the current of the SET 24 is 
monitored to determine whether a change in state occurs or not. If the system had 
initially fallen into the right 18 localised state 53, it will eventually relax to the left 17 
localised state 54. Since the state was actually in a localised state prior to this transition, 

20 the change in SET current will be detected. Based on whether or not a change of state 
is detected in the SET current, the initial state of the system can be inferred. 

Readout 

25 For definiteness we consider readout of localized qubit states with the S-gate 23 non- 
zero. A general state at time / will be of the form: 

|^(0)=C^(0|^)+Q(Ol/?>. 

30 We wish to collapse the superposition to \L) or \R), vAth probabilities IC^COI^ and 
I (Op respectively. By turning the SET 24 "on" the P-P"^ system will decohere 
strongly due to the fluctuating source-drain current and make a transition in time 7set 
« Tdissipaiion to a Statistical mixture of the localized eigenstates. In general this initial 
change in charge density will not be detected in the SET because one of the amplitudes 

35 can be arbitrarily small. Leaving the SET "on" for a further period of order of Tdisstpation 
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will allow the system to relax to the \L) state if it was in the \R) state after time Tset. 
Since the state was actually in a localised state prior to this transition the change in SET 
current will be detectable. 

5 Pakes et al have calculated SET sensitivity to a single electron moving over 20nra in 
the vicinity of the SET giving rise to a 0.36% change in the single electron charging 
energy. 

10 Initialization 

In the readout scheme as shown in Figure 16 the S-gate is set such that the |o) = \L) is 

the lower state into which the dissipative process moves the system. The process of 
initialization to the |0) state is then essentially the same as the readout process. By 

1 5 monitoring the SET current it would be possible to confirm that initialization had been 
successfully completed. 

A similar process of pre-initialization of the charge qubit itself after imbrication of the 
two donor system would be required to remove one of the electrons. The S-gate 
20 voltage would have to be higher in order to ensure an ionization event to the continuum 
takes place. 

Nanofabrication of Donor Arrays 

25 While the use of phosphorus as an electron donor within a silicon substrate is standard 
in the microelectronics industry, controlled doping of a semiconductor device at the 
atomic level, as required for the charge qubit described here, is yet to be demonstrated, 
although recent developments in the fabrication of Si:P spin qubit devices indicate that 
such precision doping is close at hand. As with Si:P charge qubits, Si:P spin qubits 

30 require buried arrays of phosphorus donors with inter-donor spacings in the range 20- 
50nm. Due to the weaker nature of the electron-mediated nuclear spin interaction, 
nuclear spin devices require very small donor spacings of order 20nm .and may also 
require positioning accuracy of up to one atomic diameter (O.Snm) In contrast, a P-P 
charge qubit with an iw/m-qubit donor spacing R=60nm still has a single-qubit 

35 operation time of IMHz and this can be made faster with the application of a positive 
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B-gate voltage. Furthennore» due to the long-range nature of the Coulomb interaction, 
the /w/cr-qubit donor spacing can be made large, possibly up to dQ-lOOnm. 

Fabrication of metallic and polymer resist features on the level 20-50nm is achieveable 
5 using existing semiconductor lithographic approaches such as electron beam 
lithography and, more recently, extra-ultra-violet (EUV) lithography. Such techniques, 
when combined with ion implantation, provide a top-down approach to array 
fabrication. 

10 The primary limitation of top-down array febrication is inaccuracy of P donor 
placement due to the statistical phenomenon of ion straggle. While such errors can 
probably be corrected by appropriate calibration of the control gate voltages a preferred 
option is to use a boUom-up fabrication approach, which offers atomic-level positional 
accuracy. Both the top-down and bottom-up approaches were developed for the 

15 fabrication of Si:P spin qubit devices and a review is given in A.S. Dzurak et al., 
Quantum Information and Computation 1, 82 (2001), 

Although the discussion below refers to controlled doping of Si with P, it is applicable, 
with minor modifications, to the doping of Si with any atom capable of acting as a 
20 single electron donor, such as other group-V elements including arsenic and 
antiminum. It is also possible to envisage the use of an alternate charge neutral 
background matrix, such as gallium arsenide, with an appropriate single electron donor, 
such as silicon. 

25 Top-Down Fabrication 

Top-down fabrication approaches involve the patterning of larger structures using 
lithographic techniques. For the creation of P-P^ donor arrays in Si we must fabricate a 
mask structure, capable of blocking the entry of P"" ions of energy l-20keV, with 

30 appropriately spaced holes allowing entry of the P* ions at the desired array positions. 
Electron-beam sensitive organic resists (such as poly-methyl-methacrylate - PMMA) 
with thicknesses in the range 20-200nm are capable of blocking keV ions. These 
organic resists may also be patterned using high-resolution electron beam lithography 
(EBL) to create via holes with diameters below lOnra. Metallic surface gate structures 

35 of thickness 20-50nm are also capable of blodcing keV ions and can be combined vdth 
organic resists in composite mask structures. Such approaches are capable of 
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producing lOnm via holes with centre-to-centre spacings (R) down to around SOnm. 
Larger spacings of R--100nm are easily achievable. It is therefore possible to create 
masks for the implantation of any desired P-P* charge qubit array. 

5 Once a mask structure with an appropriate array of small apertures has been created, 
one ion must be delivered to each cell in the array. It is therefore necessary to localise 
(or focus) the ion beam to a diameter smaller than the required donor spacing R. The 
localised beam would be moved from one cell to the next after a single P ion had been 
delivered, in a step-and-repeai approach [8]. 

10 

One way to achieve such ion beam localisation is through the use of a phosphorus 
focused ion beam (FIB) system. The same FIB system used to deliver the P ions could 
also be used to image registration markers on the chip, so as to ensure that ions are 
delivered accurately at the array sites. Alternatively the P"*" beam may be localised by 
15 collimating it with a nanomachined aperture in a mask that can be moved accurately 
with respect to the sample. If the nanomachined aperture is machined from the 
cantilever of a scanned-probe atomic force microscope (AFM) system, the AFM could 
be used to locate the array sites and then offset appropriately to position its aperture at 
the relevant cell. 

20 

Once the ion beam has been positioned at the appropriate cell, using either a FIB or 
nanomachined cantilever, the beam flux must be controlled so that only one ion 
enters the substrate. Such single ion implantation control can be achieved by utilising 
on-^hip ion impact cietectors^ capable of electrically roistering the entry of z single ton 

25 into the Si substrate. When an ion passes through an aperture at an array site and enters 
the substrate, it will dissipate its energy partly through the creation of electron-hole 
pairs in the i-Si. The electrons and holes will drift in the applied electric field, resulting 
in a current pulse in a circuit connected to the detector electrodes. Once one of these 
pulses is detected, associated circuitry can be used to blank the ion beam, limiting the 

30 implantation to just one ion. P^ ions of 10-20keV, the energy range reqmred to implant 
donors at depths of -lOnm, will produce somewhat smaller pulses, due the increased 
nuclear stopping coefficient for heavier ions. . 

In order to create the ion-stopping mask containing a nanoscale aperture array, a variety 
35 of approaches may be employed. We will discuss these with reference to the 
construction of a very simple array, containing just two donors (ie. a single charge 




32 

qubit). Here we consider three approaches, with an increasing degree of accuracy and 
corresponding process complexity. 

5 An example of a one qubit, two P donor device, febricated by the above method, is 
shown in Figure 21. The device has two SETs. By cross-correlating the current from 
both SETs it is possible to reject spurious signals and confirm true electron transfer 
between donors. The five gates at bottom are used as B and S-gates and the gate-width 
of 20nm allows for an inter-donor spacing, R=60nm. The alignment between control 

10 gates and SETs in this device is better than 20nm in both x and y directions. This 
device has about 100 P donors implanted at each site, however, this number can be 
reduced to one donor per site. The device demonstrates that the architecture described 
here is physically realisable. 

1 5 Using the first i^)proach, alignment of the S- & B-gates, and of the SETs, to the donor 
positions is achieved by roistering each stage of EBL to a series of alignment markers 
deposited on the sample at the start of the process. The accuracy of EBL alignment is 
system specific and depends on the electron beam spot size and type of alignment 
marker. In practical terms accuracy better than 20nm is challenging, however, there is 

20 no fundamental reason why accuracies as small as the spot size (below 2nm on some 
systems) cannot be achieved. The fact that this approach requires two post-alignment 
steps means that it will have inherent misalignment accuracies. This will mean that the 
B- and S-gate voltages for charge qubit operation will need to be detOTnined for each 
device, using an appropriate calibration procedure. 

25 

The first approach is to form the aperture array (in this case comprising only two 
apertures) in a single layer of organic resist, such as PMMA, using electron beam 
lithography (EBL). On-chip detectors are used to register the entry of just two ISkeV 
P"*" ions into the substrate. Following ion implantation the sample is annealed to a 

30 temperature of 900C or above, to remove lattice damage and activate the two P donors. 
Next, an array of metallic electrodes are fabricated above the donor positions, using 
EBL and standard metal lift-off. These electrodes serve as the S-gates and B-gates for 
the charge qubit. Any metal combination capable of lift-off may be used for the gates, 
such as a Ti/Au bilayer of thicknesses 20nm/40nm. Finally, two single electron 

35 transistor (SET) devices are fabricated close to the donor sites. These SETs will be 
used to detect the electron position. The SETs are made from Al, with AI2O3 tunnel 
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barriers, and are febricated using EBL with a standard bilayer resist 
(PMMA/Copolymer) and double-angle shadow evaporation process. 

Increased donor positioning accuracy can be achieved using a second, self-alignment 

5 approach. This approach uses the' metal S- and B-gates themselves as part of the ion 
implant mask. Firstly, on-chip ion detectors are fabricated for implant r^istration. 
Following this the S- and B-gates are deposited, using EBL and metal Uft-off. The 
sample is then coated with a layer of organic resist (such as PMMA) and a single cross- 
line is exposed in the PMMA using EBL. Following development of the cross-line, 

10 two apertures are produced for ion implantation. Once again, the on-chip d^ectors are 
used to limit, implantation to only two ions, and then the sample must be annealed to 
remove damage. Because this anneal step occurs with the control gates in place it 
places a stringent requirement on the choice of gate-metal. The metal must have a very 
low diffusion rate through the thin (<5nm) SiOz barrier layer at high temperatures of 

15 order 900C. This issue is also of critical importance for the fabrication of MOSFETs 
with very thin gate oxides for VLSI applications and a^ number of candidate metals, 
including W and Ta, have been proposed which have low diffiisivities at high 
temperature. FoUowing the anneal step two SETs are febricated on tiie device using an 
identical process to that described for {q>proach I. Here the SETs may suffer some 

20 misalignment, but this will only modify their charge detection senativities. 

In a third approach, perfect self-alignment between implanted donors, control gates and 
readout SETs may be achieved using a trilayer resist process, described previously in 
T.M. Buehler et al., submitted to Nanotechnology. Once agwn. tiie first step involves 

25 the &brication of on-chip ion detectors. Next, a single line is exposed using EBL in a 
single layer non-organic resist, such as PMGI, and tiien developed. This line defines 
the locus of points where the implanted donors will lie. Following this, a standard 
bilayer resist (PMMAyCopolymer) is coated above the patterned PMGI, and a second 
pattern for the gates and SETs is exposed using EBL. The bilayer is tiien partially 

30 developed to open apertures for ion implantation only where the two patterns intersect. 
Follovwng controlled implantation, limited to only two donors using the on-chip 
detectors, the bilayer is fiilly developed to open a cavity in the middle layer and the 
lower layer of PMGI is fiilly removed. The gates and SETs are then deposited by 
evaporating metal through the ronaining resist structure fix)m tiiree-separate angles, 

35 using a shadow evaporation approach. In this way, the SETs are perfectly aligned to 
the control gates, and hence the donors. The most established way of fabricating SETs 
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by shadow evaporation is using AI metal, and oxidizing this in a controlled way to form 
AI2O3 tunnel barriers. A drawback of this fully aligned approach is that the annealing 
step, required to remove implant damage, must be performed following the gate and 
SET deposition. High temperature anneals above 500C are likely to result in damage 
5 of thin AI2O3 tunnel barriers, while anneals up to 900C will lead to difRjsion of Al 
through the SiOj barrier layer. The development of SETs using different metals, with 
much lower diffusion would therefore be required. Alternatively, a lower temperature 
annealing process would need to be established. 

10 We note that the febrication of a two donor device using the first or second approach is 
already possible. Although the localisation of the implanting ion beam to single sites 
has not yet been achieved, the use of a standard blanket (unfocused) implant will lead 
to a 50% yield of ideally-configured devices, with one P donor at each of the two sites. 
This yield is acceptable for demonstration purposes and makes possible immediately 

1 5 the demonstration of a single charge qubit. 

Bottom-Up Fabrication 

The bottom-up strategy utilises both scanning tunnelling microscopy (STM) and 
20 molecular beam epitaxy (MBE) to create a P qubit array in Si, and does not suffer firom 
the problems of ion implantation straggle and substrate damage that are inherent to the 
top-down approach. The bottom-up strategy is a four-stage process. Firstly, a clean 
Si(00l) surface is passivated with a monolayer of H. This "H-resisf layer is 
subsequently patterned at the level of removing individual H atoms from the surface by 
25 applying controlled voltage pulses to an STM tip in tunnelling contact with the surfece. 
The third stage of the process is to expose the pattwned H-reast surface to phosphine 
(PH3) gas, allowing individual PH3 molecules to bond with the exposed sites of bare Si. 
The final step in the process is to use MBE to overgrow layers of Si and encapsulate the 
febricated P array. Such STM based hydrogen lithography has previously been 
30 proposed for the assembly of atomically-ordered device structures, in particular for the 
fabrication of single-electron integrated circuits. 

A completely defect-free area of a low defect density Si(001) surfiice is optimally- 
prepared using a standard thermal treatment, described previously. In the Si(001)2xl 
35 surfece there are rows of paired Si-Si atoms (dimers). After exposure to atomic H at a 
substrate temperature of -600 K a near uniform coverage of the monohydride phase. 
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where one H atom bonds to each Si atom, occurs. Should a combination of 
monohydride and dihydride phases occurs (i.e., one monohydride dimer adjacent to a 
single dihydride species) we get a H:3xl reconstruction. A surface of purely 
monohydride phase can be formed by using a slightly higher substrate temperature 
5 during H-dosing. 

An area of a hydrogen passivated Si(OOl) surface can be patterned using an STM^tip. 
In particular, two hydrogen atoms can be removed from this surfece by "pulsing" an 
STM tip to -7 V and 1 nA for ~1 ras at two locations above the H passivated surfece. 
10 These pulses have caused the desorption of two individual H atoms from the surfece, 
exposing two surfece Si atoms. Each of these Si atoms has an electron in a "dangling 
bond" oibital. The dangling bonds present reactive sites for the adsorption of chemical 
species such as phosphine - in fact the sticking coeflficient of PH3 to the clean silicon 
surface is very close to unity. The effectiveness of the H resist as a barrier to phosphine 
15 adsorption is demonstrated by the uniform H coverage after PH3 dosing except at the 
previously desorbed sites. There is a characteristic increase of -0.05 nm in the height 
of the protrusion after phosphine dosing. This increase, which was seen reproducibly at 
a number of sites, confirms the adsorption of a PH3 molecule and corresponds to the 
difference between the exposed silicon dangling bond and the adsorbed phosphine. 

20 . . 

The final step, currently underway, is the encapsulation of the phosphorus qubits m a 
crystalline lattice of silicon. The main difficulty in this step is to ensure the P atoms 
incorporate into the Si crystal and remain in an ordered array. The most direct route to 
achieving this is to desorb the H resist from the surface by heating to -825 K fijUowed 

25 by epitaxial Si growth over the P array. A possible concern with this level of heating is 
the potential to induce lateral surfiice difRision of the P atoms in the array. However 
there is evidence that the phosphorus atom in PHc (x=2,3) molecules, which has a 
single bond to the silicon, incorporates into the silicon surfiice over the lower 
temperature range of 620-720 K. 

30 

Fx perimental demonstration of single qubit rotation 

We now describe an experimental demonstration of the ability to perform and measure 
using SETS, a single qubit rotation with a P-P* charge qubiL This experiment is also 
used to detennine the coherence time T2 of the qubit. It is a direct analogue of the 
35 landmaric experiment perfisrmed by Nakamura et al [Y. Nakamura, Y.A. Paskin and 
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J.S. Tsai, Nature 22S. 786 (1999)] to measure quantum oscillations in superconducting 
qubits. 

The experiment requires an ability to measure the transfer of a single electron charge 
between two P donors on a timescale shorter than the relaxation time of the system, 
5 Tdi«ipa«ion=T,. This measurement can be made significantly more reliable by using the 
cross-correlated outputs of two SETs. as described by Buehler et al.[T.M. Buehler et 
al., cond-mat/0207597 (2002)]. An example of a twin-SET device used for such cross- 
correlation measurements is shown in Figure 18. Figure 19 shows data obtained using 
this device, which shows that both SETs respond to the transfer of single electrons 
10 between two metallic islands located between the SETs. Of particular note is the fact 
that the SETs are measuring the charge transfer in a "single shot" measurement, and 
that the measurements are occurring on the microsecond timescale. 
The limitation on the speed of SET charge position read-out is set by the shot-noise 
within the SET, and corresponds to a minimum required read-out time tmn, given by 
1 5 (/ min)"^ = 1 .2xl0-*c/ A^. where is the change of chaige that is detected. Based on 
measurements obtained by Buehler et al (2002), we find ftnin=10ns for our P-P* 
system. This is confirmed experimentally, as shown in Figure 20, which shows that 
equivalent values of to our P-P* system, when simulated using the device of Figure 
18, can be detected on timescales of around 100ns. We note that we expect relaxation 
20 times Tdisripa«io« to exceed this timescale by an order of magnitude, or greater. 

The proposed experiment to perform coherent qubit rotations, and read-out the result, is 
described schematically in Figure 17. 

In this proposed experiment a single P-P"" qubit is carefiiUy controUed by the externa! S 
25 and B gates, and measured at precisely defined intervals to determine the effect of 
decoher«nce on the quantum evolution of the system. In order to describe the physics of 
the experiment we will concentrate on an ideally fabricated P-P* qubit as shown in 
Figure 3a. 

30 The experiment is carried out in several stages. 
Stage J: localisation/irtitiaUzation 

At initial time T = 0, the S-gate is set at or near the point So = - 0.3 V, which will 
35 localiM the electron in the left eigenstate |*P(7' = 0)) =|L) . 
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Stc^e 2: S-gate pulse 

If we now pulse the S-gate swiftly to zero, with a fall-time of Tm (assumed to be 
short), the system is forced into a superposition of the So = 0 eigenstates (symmetric 
5 and anti-symmetric states): 

\nT,,oHL)=-^\s)^\A)). 
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Stage 3: Rati evolution 

While So = 0 the system is allowed to evolve through its natural Rabi transitions 
1 S) <^ 1 4 for a time Tp„ise • The state becomes: 

Stage 4: S-gate pulse 

The S-gate is pulsed swiftly back to the original T = 0 value over a time Trf». and the 
state, in terms of the localized So 5* 0 dgenstates, becomes: 



As long as Tpuise » Twuhs. we obtain essentially a state of the form: 

25 • |H'(7kl+?'pnlse + 7W«,)) = CL(rpuisc)lL)+CH(rp„,sc)l/?) • 

Stage 3: measurement 

Switching on the SET and measuring the state will give a left or tight localized state as 
3d determined by the distributions: 
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Repeating the process many times for each Tp„isc will result in a measurement of the 
5 Rabi oscillations, and the decoherence time. 

The hierarchy of time scales is important to consider. The controlling scale is of course 
the decoherence time. Although we do not know this before the measurement, we 
expect this to be of the order of 1-lOns based, on the GaAs measurement Taking this 
10 range as our lower bound, we therefore require Tjuu < Ins in order to see Rabi 
osciUations. If we use the Nakamura experimental parameters as a guide, we can set 
TRaM =100 ps and use fest pulse generators with Trfse/feii = 30 ps. 

To the extent that we have described single qubit dynamics in section 2 a time 
15 dependent analysis can be carried out on the results expected from realistic pulse 
shapes. Furthermore, the first P-P* devices will not be ideally or symmetrically gated 
and will require detailed characterization, possibly using an extra S-gate. 



It will be appreciated by persons skilled in the art that numerous 
variations and/or modifications may be made to the invention as shown in the 
specific embodiments without departing from the spirit or scope of the invention 
as broadly described. The present embodiments are, therefore, to be 
considered in all respects as illustrative and not restrictive. 

Dated this nineteenth day of August 2002 

Unisearch Limited 

Patent Attorneys for the Applicant: 

F B RICE & CO 
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